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CN ■ Abstract 

Oh> The Cauchy problem is considered for the scalar wave equation in the 



Schwarzschild geometry. Using an integral spectral representation we de- 
rive the exact decay rate for solutions of the Cauchy problem with spher- 
ical symmetric initial data, which is smooth and compactly supported 
CNj ' outside the event horizon. 



q-i 1 Introduction 



In this paper we study the decay rate for solutions of the scalar wave equation 
with spherical symmetric initial data in the Schwarzschild geometry, which is 
smooth and compactly supported outside the event horizon. We prove that 
these solutions decay at the rate t -3 and t~ A for momentarily stationary initial 
£f-^ • data, respectively, as it was earlier predicted by Price 16J, though not rigorously 

proved. In |15j we have already shown pointwise decay for solutions of the same 
kind of initial data not necessarily spherical symmetric. To this end, we have 
derived an integral spectral representation for the solutions applying Hilbcrt 
0^ ' space methods in terms of special solutions of the Schrodinger equation the 

so-called Jost solutions. 

In order to set up some notation, recall that in Schwarzschild coordinates 
(i, r, if), the Schwarzschild metric takes the form 



ds 2 = gij dx l dx^ 



-l 



= (l-— ) dt 2 - (l- — ) dr 2 -r\de+ wtQdq?) (1.1) 



r J \ r 

with r>0, 0<i9<7r, < f < 2ir. The metric has two singularities at r = 
and r — 2M . The latter is called the event horizon and can be resolved by a 
coordinate transformation. We consider the scalar wave equation in the region 
r > 2M outside the event horizon, which is given by 



-g dx l \ dxi 



U4> := gVVtfrf = —— ( _ ) = (1.2) 



where g denotes the determinant of the metric gij. We now state our main 
result. 
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Theorem 1. Consider the Cauchy problem of the scalar wave equation in the 
Schwarzschild geometry 

□0 = 0, {<j> Q ,id t <j> Q )%r,x) = $ (r,x) 

for smooth spherical symmetric initial data $o G Cfi°((2M, oo) x S 2 ) 2 which 
is compactly supported outside the event horizon. Let $(£) = (<fr(t) , idt(j>(t)) € 
C°°(R x (2M, oo) x 5 2 ) 2 be the unique global solution which is compactly sup- 
ported for all times t. Then for fixed r there is a constant c = c(r, $o) such that 
for large t 

\m < £ ■ 

Moreover, if we have initially momentarily static initial data, i.e. dt4>o = 0, the 
solution <p(t) satisfies 

\m < £ • 

There has been significant work in the study of linear hyperbolic equations 
in black hole spacetimes. The first major contribution in this topic was made 
in 1957, when Regge and Wheeler studied the linearized equations for pertur- 
bations of the Schwarzschild metric [T7] . This work was continued in [521 US] > 
while more recently the decay of the perturbation and all of its derivatives was 
shown in [T^] using a theorem by Wilcox. By heuristic arguments, in 1972 Price 
[IB] got evidence for polynomial decay of solutions of the scalar wave equation 
in Schwarzschild, where the power depends explicitely on the angular mode. In 
1973, Teukolsky [T<|] could derive by means of the Newman Penrose formalism 
one single master equation that describes in the Kerr background the evolution 
of a test scalar field (s = 0), a test neutrino field (s = ±1/2), a test electro- 
magnetic field (s = ±1) and linearized gravitational waves (s = ±2). Here, 
the parameter s is also called the spin weight of the field. Note that it is a 
quite complicated task in the case s ^ to recover all the components of the 
corresponding field from a solution of this equation. For further details see 
[3l [24] . In two subsequent papers [20l [21] , Teukolsky and Press discussed the 
physical consequences of these perturbations. Note that the rigorous analysis of 
the equation remains a quite subtle point, though any linearized perturbation is 
given by this equation. For instance, in the case s ^ complex coefficients are 
involved, which makes the analysis very complicated. Hence, until now there 
are just a few rigorous results in this case. In [7] local decay was proven for 
the Dirac equation (s = h) in the Kerr geometry (in the massless and massive 
case). Moreover, a precise decay rate has been specified in the massive case 
[5] . More recently, there has been a linear stability result for the Schwarzschild 
geometry under electromagnetic and gravitational perturbations [S]. This re- 
sult relies on the mode analysis, which has been carried out in [23j . More work 
has been done on the case s = 0, where the Teukolsky equation reduces to the 
scalar wave equation. In the Schwarzschild case, Kay and Wald [2] proved 
a time independent L°°-bound for solutions of the Klein-Gordon equation. In 
[5], a mathematical proof is given for the decay rate of solutions with spher- 
ical symmetric initial data, which has been predicted by Price [16j . which is 
not sharp, however. For general initial data, the same authors derived another 
decay result [Bj. Pointwise decay in the Kerr geometry was proven rigorously 
[TU1[TI]. Furthermore, Morawetz and Strichartz-type estimates for a massless 
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scalar field without charge in a Reissner Nordstr0m background with naked sin- 
gularity are developed in [18] . And in [2] a Morawetz-type inequality was proven 
for the semi-linear wave equation in Schwarzschild, which is also supposed to 
yield decay rates. 

In this paper we first recapitulate the framework and some notations of the 
foregoing paper [13]. Afterwards, we give an explicit expansion of the Jost 
solutions <f) of the Schrodinger equation which also were derived in [IS]. At the 
end we show how to derive the exact decay rate out of this expansion. 



2 Preliminaries 

We usually replace the Schwarzschild radius r by the Regge- Wheeler coordinate 
ael given by 

u(r):=r + 2Mlog(^-l) . (2.1) 

After having separated the angular modes I, m using spherical harmonics, it is 
convenient to write the Cauchy problem in Hamiltonian formalism 

*ft* = H*, *| (=0 = *o (2.2) 

where = (ip, ii9tV') T is a two component vector representing the wave function 
and its first time derivative and H is the Hamiltonian 



1 
+ 



(2.3) 



with the potential 



^) = { 1 -—){— + ± ^r 2 ) ■ ( 2 - 4 ) 



Constructing the resolvent of the operator H and using Stone's formula we 
have derived an integral spectral representation for the solutions of the Cauchy 
problem of the following form 

-- I e+* (f iJ ( \ M Mv)dv) du> , (2.5) 

where the integrand is in L 1 with respect to u>. 

At this point, the functions 0, <p play an important role. These functions are 
a fundamental system of the Schrodinger equation 

(-% + V u (u))<f>(u) = (2.6) 

with the potential 

9 / \ of 2M\ f2M 1(1 + 1)\ , , 

V u (u) = -J 1 + Vi(u) = -uj 2 + h - - ( — + ^r 1 ) (2.7) 

with boundary conditions 

'^(«) = 1, lim (e-^ u 4,(u)Y = (2.8) 
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lim e iuu 4>Ju) = l, lim (e iuu ^ w (u))' = . (2.9) 



in the case Im (w) < 0. We derived these solutions using the corresponding 
integral equation, the so-called Jost equation, which is given by 

<j> u (u) = e*" u + f -wn(u(u-v))Vi(v)<f> u (v)dv, (2.10) 

in the case of boundary conditions at — oo (an analog equation is considered for 
the boundary conditions at oo). Now, we have constructed the solution <j> with 
the series ansatz 

oo 

&, = $><, fc) , (2.11) 

together with the iteration scheme 



/U i 
-sm(w(u-«))V5(t#|f>(t;) dv 
-oo w ) 



(2.12) 



Using this we have proven that the solutions <fru(u) are analytic with respect 
to uu for fixed u in the region Im (w) < 0. Moreover <j> u can be analytically 
extended to the region Im(w) < jjg, whereas the solution u l <f> u can be only 
continuously extended to the real axis. Thus, in order to obtain the exact decay 
rates it is important to analyze the behavior of W with respect to w at the real 
axis in more details. 



3 Expansion of the Jost solutions W 

Since the w-dependence of the Jost solutions (j> u plays an essential role in the 
analysis of the integral representation, we show in this section a method to 
expand these solutions at the critical point u) = 0. We start with an explicit 
calculation: 

Lemma 3.1. For all u > 0, uj e R \ {0}, e > 0, q e N and peN, 



/>OG 

J u 



m— 



f_1 f_1 \m+l m 

( J +1 log m+1 [(2iw + + J2 log* [(2i W + e)u] 



fe=0 

oo 



w/iere i/ie coefficients Ck involve the coefficients do, ...,a q of the series expansion 
of the T-function at 1 — p. 
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Proof. In order to prove this, we write the integral as A-derivatives, 



e -2iu, a -e >g g (*) dx = ^_ F (X) 

6 xP X ~d\i p[A) 



(3.2) 



A=0 



with the generating functional 
F P (X) -- 



OC 1 

(-2iu,- e )x_!_ dx = U -P+X+1 



x" 



-A 



= ( — 2iu> — e)i 



V P- 



dv 



where in the last step we introduced the new integration variable v = — . In the 
following we will write z = (2iu) + e)u for reasons of convenience. The integral 
on the right hand side is also known as the Exponential Integral E p -\(z) with 
the series expansion 



E p . x (z)=T(l-p + \) zP-^ 1 -^ 



(-l) fe 



(k-p + \ + l)k\ 



for small A ^ [as a reference cf. Using the series expansion of the T- 

function at 1 — p e Z \ N, where the T-function has a pole of first-order, we 
obtain 



F P (X) = u 



(-I)"" 1 
(p-l)!A 



^ a„A™ 



ri=0 



E 

fe=0 



(-1)* 



(fc-p + A + l)fc! 



-p+A+1 



-i [ (-IF" 1 



(p-iy. 



+ Z~ 



E (fc - + A 4- 1 )/■'! 



Using z — e gz , we immediately get the formulas 



(3.3) 



dF 1 



z~ x - 1 



d\ r 
d m 
d\ m 



r(" A ) 



A=0 



A=0 



A=0 



(-l)" +1 log" +1 (z) 
n+1 

log m («) 
(-l)" l log m (z) 



one directly verifies the claim setting (|3.3|) in 



□ 



Directly in the same way, one proves an analogue lemma for the case p G 
Z\N: 

Lemma 3.2. For all u > 0, uj 6i\ {0}, e > 0, q e N and p e Z \ N, 

.A, 
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' q 



E (m) loga ~ f »i (2^ + er 1 E c ^ m ) 1 °g fc [( 2 ^ + £ H 

fe=0 

„ HX . (-l) fe (-l) m m! 



i m 

m=0 v 7 k fe=0 



k=0 



l Y, ,u ,L?iu [(a^ + eH* (3-4) 



where the coefficients Ck involve the coefficients clq, a q of the series expansion 
of the T -function at 1 — p. 

Compared to Lemma l3.1| here the logarithmic term is of lower order due to 
the fact that the Gamma-function has no singularity for positive integers. 

In order to apply this lemma to our integral representation, we have to derive 
an asymptotic expansion for the potential Vi(u) at +oo. Therefore, we have the 
following 

2M\ ( 2M 1(1 + 1)* 



Lemma 3.3. For the potential VAu) =11— , , 

\ r(u) J \r(uy r(u) 2 

have the asymptotic expansion 

*H") = 2^2^ c pi— -p— + c k+iM-i k+1 +o\ uk+1 , (3.5) 

p=2 q=0 \ J 

as u — > oo , with k > 2 and real coefficients c pq , where e.g. the first coefficients 
are given by 

c 20 = 1(1 + 1) , C 3 1 = 41(1 + 1)M , 

c 30 = 2M - 2Ml(l + 1)(1 + 2 log(2)) - 4Ml(l + 1) log(M) 
c 42 = 12Z(Z + 1)M 2 , 

C4i = -UI 2 (-3 + l(l + l)(5 + 81og(8)) + 61(1 + 1) log(M)) , ... . 

Furthermore, in the case I = the coefficients c„.„_2 vanish. 

Proof. First we have to find an expression for r in terms of the Regge- Wheeler 
coordinate u. Remember that u = r + 2Mlog(2^y — 1), which is equivalent to 

f 2M~ 1 = | — 1 ] p~TM~^ 

\2M J 

In order to resolve this equation with respect to r, we use the principal branch 
of the Lambert W function denoted by W(z). This is just the inverse function 
of f(x) — xe x on the positive real axis. [As a reference cf. [1].] Hence, we obtain 

r = 2M + 2M W(e-™- x ) . (3.6) 

Moreover, for W we have the asymptotic expansion 

oo oo 

W(z) = log z - log(log z) + E E Ckm ( lo g( lo S z )) m+1 ( lo § z^™" 1 , (3-7) 



fc=0 m=0 



as z — > oo. Here, the coefficients Ckm are given by Ckm = "^~r(~l) fe 



k + m 
k+1 



where 



k + m 
k+1 



is a Stirling cycle number. In particular, applying this expansion 
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to (|3.6p , we get the series representation 

i.2M 



[XL f U 

^-l-logf — -1 



oo oo 



^2 Ckm ( lo gf-r^T - I 



k=0m=0 



\2M 



m+1 



f— -i 



-k—m— 1 



This allows us to expand the powers -\, and ^ to any order in u/2M — 1 
using the method of the geometric series. Together with the expansion 



u ( 2M\~\ , , , ^ 1 (2M 

— 1 = log U -log2M - V- 

"M \ u I ' n \ u 

n—l v 



which holds for u > 2M, the result follows. □ 

These two lemmas let us expand the solution (j)^ (u) in the following way. 

Lemma 3.4. For I = 0, u) G M. \ {0} and fixed u > 0, i/ie fundamental solution 
4>u,(u) can be represented as 



<i>u(u) — e + g (uj 7 u) + 2iivlog(2iuj)g 1 (uj,u) + 2iug2{u 1 u) , 
where the functions go,gi and gi are C (M) with respect to to. 

In order to prove this, we need the following lemma 
Lemma 3.5. For all u £ C and n € No, 



(3.8) 



9™ I — sin u 

u 



Moreover, i/w^0 and v > u > 0, 



1 



C[-sin(w( U - V ))] 



on+l 
< J JImul 

- i + N 



^ C(n) «" + u |i mtl ;|4-ttlmh) 



(3.9) 



(3.10) 



1 + 

/or some constant C(n), which is just depending on n. 

Proof. In the case \u\ > 1, Eulers formula for the sin- function inductively yields 



(i + H) 



9™ ( — sin u 

u 



< 2" +1 e' Imu ' 



For \u\ < 1 we rewrite (l/u)sinu as an integral, in order to obtain the 
estimate 



(i + H) 



i)'.'. ( — sin-u 
u 



= (1 + H) \ I (ir) n e mT dr 



< 2e |Imu| 



which shows the first claim. As a consequence, we get for uj ^ and all n e N 
the estimate 



d:[-s\n(uu) 

id 



UJU 



< 



I Im (uju) | 



1 + \luu\ 

In order to show (|3.10[) we use the identity 

- sm{uj(u -v)) = - (sin(uju)e tuv - sin(ujv)e^ u ) 

UJ UJ 



(3.11) 



(3.12) 



and apply (|3.9p . n = 0, 

— sin(ui(u — v)) 

UJ 

2\u\ 



< 



sm(uju)e iu)v \ + |sin(ww)e i " u | 



^ ~ i ~ i Im uj \ ^—vlmm 



2\v\ 



Im a; | ^— u Im uj 



1 + \uju\ 1 + \ujv\ 

Due to the assumption v > u > 0, we know that |t>| > \u\ and thus 

2 M 2|u| ■ ■ 

< ; r, u\ lmu>\ + ulmw > v\ Imw| + ulmw 



(3.13) 



1 + |um| ~ 1 + \ujv\ 

Using these inequalities in (|3.13p the claim follows for n = 0. 
Once again using the identity 13.121 we get 

d w ( - sm(uj(u - v))) < - (sin(Luu)(-iv)e- iuJV - sm(wu)(-m)e~ iwu ) 

\UJ J UJ ' 



1 



1 



Using the estimates p. lip and (|3.9p for n = together with the assumption 
v > u > 0, we see as before that the first term is bounded by 



1 



(sin(uu){-iv)er l " v - sm(ujv)(-iu)e- iuu ) 



^ lmuj\ J ru ImoJ 

~ 1 + \ujv\ 



For the second term we use ([3.11)1 



d u - sin( WM ) e~ tuJV - d u - sin(wu) e""" 



Aip Alp 

^ttl Im a;|+-u Im u _j_ Im u | +^ Im it; 



l + |wu| l + |wu| 

and obtain due to the assumption v > u > 

^ ^^|Imuj|+Tt Inm 

~ 1 + \wv\ 

Thus, we have shown (|3.10p for n — 1. We proceed inductively to conclude the 
proof. □ 

Note that the estimate (|3.10p remains valid in the limit ^ ui — > for all n, 
because 



lim^(-sin(w(«-t;))) = 

lj^Q \ OJ ' 



(_l)«/2^_( u _ w )n+i if n even, 



if n odd. 
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Proof of Lemma \3.4\ First, remember that the solution <f>u(u) is given by the 
perturbation series 

oo 
k=0 

where the summands follow the iteration scheme 

r°° i 

^£ 0) («)=e- fc " i , $ +1 Hu) = - / -Mu(u-v))Vo(v)<l>W(v)dv, (3.14) 



with potential Vb(u) = 1 



2M\ 2M 



r(u) J r(u) 3 



According to Lemma [3.3[ this poten- 



C30 



tial can be represented for large u as Vo(u) = —r+hOu), with h(u) = O 

u 6 

Next, we split this iteration scheme up. To this end, we define 



4£\u) := 



and analogously, 



& x) («) 



-sm(uj{u~v))h(v)e" wv dv 



- sm(u(u - v))^-e~ luJV dv 



\ogu 



(3.15) 



(3.16) 



Thus, obviously (j>jj\u) = (u) + ffl (u). Now we iterate these two functions 

f°° 1 

^ +1) {u) := - / -mi{uj{u-v))V Q {v)^\v) dv , fc > 1 , 



(3.17) 



analogously for Hence, we have the formal decomposition 

oo oo 

= e— + £ 0( fc ) («) + ~€ ] («) ■ 
fc=i fc=i 

Both series are well-defined. In order to show this, we use the bound 



— sm(cj(ii — v)) 



< 



4M 



1 + 



(3.18) 



from Lemma 13.51 for real ui [Note that this estimate is also valid for the case 
v > u > 0]. Hence, we get inductively the estimates 



fc! 



\^Hu)\<R»(u)^^, 
for all fc > 0, where the functions R, R and P are given by 



ilW := 



-it; 



1 + \u)\v 

Av 
1 + \oj\v 



C30 
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P u (u) 



4v 



\V (v)\ dv . 



L 1 + \u\v 

Thus, the series ^ <j>ffl (u) as well as ^ cj>ffl (u) converge locally uniformly with 
respect to u and u>. In the next step we show that, for fixed u > 0, 4>ffl ( u ) 1S 
C 1 (R) with respect to uj. To this end, it suffices to prove that each summand 
4>£\ k > 1, is C 1 and that the series J^c^^w converges locally uniformly in 
uj. Due to the estimates (|3.10|) . (|3.18[) . we have the inequality 



9. 



-sin(uj(u~v))h(v)e- luJV 

UJ 



< 



< 



12 v 2 



1 + \cj\v 



h(v] 



Av 2 



-h(v 



1 + \lu\v w 1 + 



16u 2 



(3-19) 



Hence, the second term is an integrable bound, uniformly in uj, for the first 
derivative of the integrand. It follows that 4>u\u) is C 1 with respect to uj, 
bounded by 



\dJW( U )\ < 

Together with the estimate 
1 



16v 2 



u i + M« 



\h(v)\dv=:R^(u) . 



4it / u 1 + u 
one shows inductively that 4>^ +1 \u) is C 1 with respect to u>, bounded by 

\dJ^\u)\<RU(u) {4P ^ ))k . 

This yields that the sum ^ c^^L converges locally uniformly in uj. Hence, 
the sum <j^\ u ) is C 1 (R) with respect to uj. According to the decomposition 
p,17p . it remains to analyze the a;— dependence of J2^\u). To this end, we 
compute the first summand: 



»(«) 



/ 1> 3 



1 



c 30 , 1 », 



2iw J u v 
Integrating the second term by parts, we obtain 

^ f — iuu iuju ^30 —2iuu i iuju 

2iuj\ 2u 2 2u 2 + 

/•OO 

2u^ 



yd 



C30 

-2u 2 



The series expansion of Lemma 13.11 in the limit e — > yields 



= ^e iwu <j2^(log(2^) + c ) 
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^ (k-l)k\ 



(2lU)U) 



(3.20) 



Intuitively, the only term which is not C 1 is the term involving 2iuj\og(2iuju). 
More precisely, defining 



^(u) := ffi\u) - C8oe iut, iwlog(2i«u) , 
and iterating this by 

f°° 1 

^ +1) {u) :=- - sm{w{u - v))V (v)^(v) dv,k>l, 



(3.21) 



we show next that the sum ^ipu (u) is C 1 with respect to u>. By definition 
this holds for the initial function -ipffl (u) . In order to prove this for the sum, we 
apply the same method as above. To this end, we need good estimates for the 
initial functions $y (it) and (u). Estimating the integral representation 

of $\u), we obtain for arbitrary u > and 



+c 30 e^iw log (2iwu) = C J («) 



< 



C30 



2w 2 



dv = 



£30 

2w 



On the other hand, looking at the series in (|3.20p , we obtain for all u < the 
estimate 



I 7(i)/ ni c 30 i WU r o . 1 v-^ (~!) fe / - \k\ 

W{u)\= —e {2 M „-- ^ ^_^ (2iWU ) } 

with a suitable constant c. Thus, we get for all u > and cj € K the estimate 

c 

where c is chosen suitably and l.(.) denotes the characteristic function. In order 
to estimate the derivative d UJ 'tpw\u), we use in the domain u > j^j, \u>\ ^ 0, the 

following bound for d w 4>u\u) [see also (|3.19[) ]. 

16v 2 



\^(u)\ < ^+cH|log(2to«i)|l [ i )0o) («) , 



< - , (3.22) 

it 



(3.23) 



W (0£ X) («) 



< 



< 



u 1 + M« 



16 



IC30I , , 16c 30 

— 5- «w < 1 1 < 16 c 30 



Together with the analogon to estimate (|3.22p in the region u < j^j , we obtain 
the bound 



\d u 4>jp («)|<c+ 5(1 + u\u\)\ log(2ium)| 1[ 1 i0o) (u) 



(3.24) 



where «>0,wel and c is an appropriate constant. For reasons of simplicity, 
we choose c = c such that both inequalities (|3.23[) . (|3.24p hold. Using these 
inequalities, we show by induction, in the same way as above, that tJj^ (u) is 
C 1 with respect to u and obeys the estimates 



cP u (u) 



k-l 



u (k — 1)! u 



r-(M)- 



(k-2) 



(3.25) 
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for all > 2, it > and w€R, where r is given by 

KM) := |Vo(«)|| log(2i^)| dv . 

Due to (|3.25[) . (|3.26[) , the sums ^ V>£, (u) and ^ i u ) converge locally uni- 

formly in to. Hence, we conclude that ^ $w (it) is well defined and continuously 
diffcrcntiablc with respect to uj. 

Thus, it remains to look at the term we get by the iteration of 

i?W(u) := c 30 e iuu iujlog(2ioju) = ic 30 wlog(2iw)e iam + ic 30 uj log(u)e luJU . 

To this end, we split up the iteration, exactly as we did for the iteration of 
4>u (u), i.e. we define 

f°° 1 

^ 2) (u) := -/ -sw(w(u-v))h(v)#W(v)dv, 

Ju w 
roo 1 

^(u) := -/ -sm(oj(u-v))^W(v)dv, 

Ju U V d 

and iterate these functions, 

$l k+1) (u) := - / -sm(oj(u-v))V (v)-d^(v) dv , fc > 2 , 

analogously for $u (u). Next, in exactly the same way as for <j)^ k \ one sees 
that 

oo 

J2^ ) (u)=2iujlog(2iu > )f 1 (cj,u) + 2icjf 2 (u;,u) , 

k=2 

where fi(.,u) and f%{.,u) are C 1 with respect to uj. Finally, by an exact calcu- 
lation 



oo 



1 



( 2 ) (u) = ic 2 30 lu log(2iw)e-^ u / e 2 ^" ^ dv 

Uc*uje-^ u I e 2 ^ + ip^l ,/V 



V 4w 2 2w 2 



«/ It 

together with the series expansion of Lemma ETT1 in the limit e — > we obtain 
^ 2) («) = ~»4> uj(1 + 21og(2i W ))e-- " 

x (-2iw)(log(-2iwu) + co) -u- 1 J! _ lU! (~2»^) A 



+ ^c 2 we— 



E log 1 - m («){(-2^)x 

m=0 ^ ' ^ 
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^^-j log m+1 (-2iwu) + ^ c fc log fc (-2^w) 

fc=0 



^ (fc-l) m+1 fc! v ' 

k=0,k^p-l y ' 

Proceeding in the same way as for V'w (w) [i.e. we omit the logw-terms in 
the square brackets, and iterate these functions], we again get terms of the form 

2iLo\og{2iio) fs{io, u) + 2iu)fi{u)^ u) 

with continuously differentiable functions fs(.,u), fi(.,u). So after simplifica- 
tions there remain terms of the form 

(2iujf \og s {2ito) log r (-2iw) log*(iOe- iam . 

These are obviously C 1 with respect to to and so is their iteration, due to the 
fact that the additional w-order yields directly integrable bounds for all to. This 
completes the proof. □ 

Note that one can apply this idea of the proof to the case I > 1. This yields 
a similar result but also requires much more complex calculations according to 
the construction of the Jost solutions <p u [cf. jT5J Section 5]]. 

4 The decay rate for spherical symmetric initial 
data 





f e~ wt 1 


/ Im( 




K \ 


J supp 



According to the integral representation 12751 the solution of the Cauchy problem 
for compactly supported smooth initial data £ Cq°(R) 2 has the pointwise 
representation 

' kitofa<v) \ f U 1 \ \ ^ 

Our goal is now to use the Fourier transform Ij2.5|> . in order to get detailed decay 
rates. To this end, we have to analyze the integral kernel, hence essentially 

lJ^i$^M) . (4.1) 

Since we already know that <f> u is analytic on a neighborhood of the real line, 
it remains to understand <p u at the point u> = 0. To this end, we want to use an 
expansion as in Lemma l3.41 The problem is that this expansion is not sufficient 
for this purpose. Thus, we apply a similar method in order to gain 

Lemma 4.1. For I = 0, to S R \ {0}, n > 3 and fixed u > 0, we get for the 
fundamental solution 4>^(u) the representation 

n 

^(u)=e- iuju + g (u),u)+ (2iu;y\o^(2tuj)log k (-2iLj)g ljk {u; 7 u) , (4.2) 

i>j+k=l 

where the functions go,gijk S C"(R) with respect to to. 
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In order to prove this, we need the following lemma. 

Lemma 4.2. Let u > 0, n € N and h S C°°(R+) be a smooth function satisfying 

f™v n+1 \h(v)\ dv < oo. 

Then: 



poo 1 

fEKu) •■=- ~ sm{cj(u - v))h(v)e- luv dv 



is C n (R) wzt/i respect to to. 
(ii) For allk>l 

roc i 

/( fe+1 )( U ) := - / - sin( W ( U - u))Vb(tO/J fc) (t;) dv 



are C n (R) wi£/i respect to lo and the series X)fc>i ^"Z^ rn < n, 
converge locally uniformly. 

In particular, ^ f^ k \u) is C ra (R) with respect to lo. 

Proof. This is shown in exactly the same way as the statement that the functions 
in the proof of Lemma [3.41 as well as the series are C 1 with respect to to. 
In order to show the differentiability up to the n-th order, we use the estimates 
of Lemma 13.51 □ 



Proof of Lemma \4-l\ Because of complex calculations we show this at first in 
the case n — 3. To this end, we split up the iteration scheme (|3.14|) of the 
fundamental solutions in the following way. According to Lemma 13. 3[ we can 
write the potential Vq as 



5 P~ 3 i_ qi 



p=3 q=Q 

where r§ is a smooth function for v > u behaving asymptotically at infinity as 
Thus > defining 

4>EHu) ■= - / -sin(w(« - v))r 6 (v)e- lM ' dv 

and for all k > 1 

~€ +1) («):=-/ - bui(w(u -v))V (v) dv , 

Ju U 

Lemma 14.21 yields that Yl $0 S C 3 (R) with respect to cj and is a contribution 
to go(u>,u) in the statement of the lemma. Thus, we have to compute the 
remaining term 



poo 1 5 p 3 

£>(u):=-/ -Bin( W (tt-«))X;E 



5 p-3 . „ 

log 9 

p=3 8=0 



C p9 < '"" dv. 
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We do this essentially in the same way as we computed the terms (jfi) , in 
the proof of Lemma [3~4l We split up the sin(a>(u — v)) with Euler's formula and 
integrate by parts and obtain 

4 p-2 „ \ 

—^ + 2^2^c pq ——\e dv , (4.3) 



-2v 2 

p=3 9=0 



where the coefficients c pq depend on the integral functions of the terms log r v /v s 
Now, we apply Lemma |3. II in the limit e \ and get 



2iu>\og(2iuiu) dk{2iuju) h 



2 I ov ' u 

K k=0 



e ^ u £4i , , Z\ og *(2iuv)+\og{2iiou)(c- £logu) 



t 



+ (c + c log it) d fc (2j 

+ e^ u i(2iuj) 3 clog s (2ium) log* 

^ s+t=l 

2 oo . 

+ ^ c log m w — ^<4(2ium)M (4.4) 

m=0 fc=0 ^ 

with appropriate constants c and dfc, which are of the form 

( — l)k( — l) m 77i! 

dfc = t; tt — rrrr • [cf. Lemma GO 

(k-p+l) m+1 k\ L 1 

Since the series-terms are obviously C 3 (]R) with respect to u), this expression of 
(jffl (u) fits into the desired expansion (|4.2p . In the next step we have to iterate 
(|4.4p . To this end, we treat each term in the curly brackets separately. We show 
this exemplarily for the first term which we denote by 

{ OO x 
2iwlog(2ium)- -^4(2ium)M (4.5) 
u fe=0 ' 



2v z 

In order to derive sufficient bounds for all u > 0, we use different methods for 
the regions \uo\u > 1 and \oj\u < 1. First, let u be such that \ui\u > 1, and by 
integrating by parts we get: 

r°° i 
a {1) (u) = e wu \ — - d 3 e- 2luJV dv 

= ^_ p -iuv J_ _ c 30 1 , 3c 30 -ium 1 / 4 r n 

4iw v? {2iuY u 3 (2iuj) 3 u 4 V ; 
12 1 



u v 5 (2iw) 3 
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Using this expression and elementary integral estimates, we get for all u > 
satisfying \u>\u > 1 the bounds 

1 



< c- 
\d u aW(u)\ < c-\- 
\d*aW(u)\ < c^- 



|S2aW(u)| < cf- 



(4.7) 



with suitable constants c. Moreover, comparing the infinite sum of (|4.5p with 
the exponential function, one directly sees that it is C 3 with respect to w. ft 
satisfies for all u > with Icolu < 1 the bounds 



d k {2iuiu) k 

I u z — ' 

fc=0 

\ U k=0 
/ 1 00 

d 3 -^d,(2zc U ) 



< 



c 
u 



c 



< cu 



(4.8) 



Using (|4.7| . ()4.8p . we verify that iterating the sum first with r$ followed by the 
full iteration with the potential Vo, we obtain a C 3 -function. For the remaining 
term cso/2e wv 2iu)log(2iuiu) we use the identity log(2iuju) = \og(2iuj) + logw to- 
gether with Lemma l4~2"l to show that the first iteration with r$ followed by the full 
iteration with the potential Vq yields a term of the form 2iuj log(2iw)/no(w, u) + 
UJ fwo( u} , u )i /liOj/ioo € C 3 (M) with respect to uj, fitting into the expansion 
Thus, it remains to compute the integral 



■ sm [uj [u — v 



p=3 (j=0 



We do this exemplarily for the term 



— sin((j(it — v))-^-a^ (v) dv 



(4.9) 



A complex calculation, using Lemma 13.11 and Lemma 14. 3[ which will be stated 
and proven afterwards, yields 

(3 (2) {u) = 2iuj\og{2iuj)e- luJu c{{~2iuj)\og{-2iuiu) - - ^ d k {-2iwu) k \ 

^ U k=Q ' 

'-^|(-2^)(-ilog 2 (-2^) 



2iioe 
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1 00 1 

+ log(— 2iu>u)(c+ logu)) - (c + clogw)- dk(-2iuju) k \ 

^ t — n J 



fc=0 

oo 



c(2iw) 2 log(-2iww) + ^^4(-2^u) fc f. , (4.10) 

with suitable constants c, dfc. Hence (3^(u) goes with (|4.2p . So far, we cannot 
finish this scheme, but if one has a close look, one sees that the most irregular 
term at u> = 0, namely 2icj log(2iw), now appears with a l/u decay, while the 
other irregularities appear with an additional w-power. Furthermore, due to the 
bounds (14. 7|) together with direct integral estimates, we obtain for all u with 
\u>\u > 1 the bounds 

v 1 1 , 1 

Of < C — r; — r 



|/? (2) H! 


< 


So,j8 (2) («)| 


< 




< 


dl^{u)\ 


< 



1 + \ui\v v 3 v 2 \oj 

1 



c^- . (4.11) 
u\uj\ 

Using in the region \oj\u < 1 for the sum-terms in (|4.10| estimates analog to 
(|4.8p . we conclude in the same way as before, that iterating these first with re 
followed by the full iteration with the potential Vq and summing up, we obtain 
C 3 — terms. We split up the remaining log-terms by log(— 2ioju) = log(— 2iuj) + 
log(u) and use Lemma 14.21 to show that applying the same procedure yields 
terms that go with (14.21) . Hence, wc have to analyze the integral 

5 p-3 

EE' 

p=3 q=Q 



/ IsinM^,))^^c P9 ^)^)(,)^, 



exemplarily we treat the term 



r°° i 

7 (3) (w):=-/ -sm(uj(u-v))^-f3^(v)dv . (4.12 



to 



Computing this expression with the same methods one sees that the term with 
2iu;\og(2iuj) decays as l/u 2 and the lu 2 log s (±2iw)-terms decay as log'(u)/u. 
With bounds analog to (|4TTT|I .(|4T8 |) the same procedure applies and yields terms 
that match with (|4.2p . Once again it remains to analyze 



isinM M -,))^^ Cp9 i^M 7 ( 3 )(,)^ 
p=3 q=a ' 



and exemplarily 

— sin(w(it — »))-tt' 3 '(") dv 



Calculating this, one checks that the 2iuj log(2iw)-term decays as l/u 3 , the 
lo 2 log s (±2iw) -terms decay as log'(u)/u 2 and the u) 3 log m (±2iw)-terms decay as 
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log n (it)/w. Applying this scheme two times more, all terms which are not C 3 
with respect to ui decay at least as log s (M)/u 3 . Subtracting these terms from 
the full term, we obtain a C 3 -term which is decaying at least as log s (ii)/zi 3 , 
according to estimates analog to (14.11|) . (|4.8|) and estimating |w| by 1/u in the 
region \oj\u < 1. So Lemma l4~2l applies for the full iteration with the potential Vq 
and we get a C 3 -term. Due to their decay, we are able to iterate the subtracted 
log w-terms also with the full potential Vq and get terms that match with (|4.2p . 
Thus, the scheme can be stopped after finitely many calculations and the lemma 
is proven for n — 3. For n > 4 we split the potential in the way 

n+2p-3 . q , k 

^) = ££c pg ^ + r„ +3 ( V ), 



p=3 q=0 



and proceed with the same calculations. In (|4.6p we have to integrate by parts 
up to the n-th order, in order to obtain as analogon to estimate (|4.7[) 



|3> (1) (u)| < A ? 



.2-1 



m < n . 



The next difference appears in the estimates (|4.11[) . These cannot be done for 
?i > 4 by simple integral estimates as a matter of convergence. Thus, we have 
to subtract from the result of the analog calculation to (|4.6|) for a^\u) the first 
n — 3 exact terms of the form 



and get for m < n 



< 



< 



LU V 6 \ ) 

— sin(oj(u — v))-^-p^(u) 



where for the first integral this can be done by elementary integral estimates, 
and for the second integral we have to integrate the subtracted terms by parts, 
as we did to obtain the estimates for (u). Keeping these differences in mind, 
we can conclude exactly in the same way as for n — 3, which yields the claim 
for arbitrary n. 

□ 

We now state the missing lemma. 

Lemma 4.3. Let u > and \ {0}. For the calculation of the iteration of 

the infinite sums that appear in the integration in Lemma \S.l\ with an arbitrary 
part of the potential, \og q u/u p , cf. Lemma \3.3\ we obtain the identity 



1 log 9 v log 5 v ±iuJV ^ 

-sm{u(u-v)) —e* > 

UJ v p v * — ' 

fe=0 



d k (±2iwv) k dv 
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-. s+q oo 

iP —2 E cl °g'( u ) E <MT2i^) & (4.13) 



s+q m+1 



+ {2iujf +t ' 2 e Tluu E dog r (T2^u) \og s+q - m (u) , (4.14) 

m— r— 1 

/or suitable constants dki,c. 

Proof. Let us denote m = g + s > and n = p + i > 4. In order to compute the 
integral on the left hand side in the lemma, we insert a convergence generating 
factor 

- / - S m(uj(u-v))-^-^e ±luJV ^ d k (±2iujv) k dv 

Ju " «" fo 

poo i i rn oo 

= limy e~ £V -sin{u{v-u))^^e ±lulv ^d k {±2iuov) k dv. (4.15) 



k=0 



In the next step we interchange the integral and the infinite sum. This can be 
done for any e > 2\lv\ by a dominating convergence argument, if one estimates 
the modulus of the sum very roughly by exp(2|w|u). Thus, the two expressions 

- sin(w(w - u))e- £V -^-e ±luv ^ d k {±2iujv) k dv 
- W ' k=o 

and 

oo poo -t i m 

d*(±2iw)* j - sin(cj(v - u))e- £V -^-e ±MV v k dv 

k=0 Ju ^ v 

coincide for any e > 2\u>\. Moreover, both expressions are analytic in e for 
Ree > 0. So by the identity theorem for analytic functions both expressions 
coincide for any e > 0. So (|4.15|) is equal to 



00 y>M j 

lim \^ d k (±2iuj) k / — sin(oj(u — u))t 
k=o Ju 



_ ev la^v e±iuv dv 



Once again we rewrite sin(w(i> — u)) with Eulers formula and integrate by parts 
[note that one has to be careful with the e-terms that are generated by this 
integration by parts, but in the limit e \ they vanish] to obtain 



oc ,.oc jn 

lim d k {±2iuj) k e TtuJU / e ±2luJV - £V ^ , 



fog' V , 

fc=0 " u i=0 



with suitable constants c arising from the integral function of log™ v/v n k . Now 
we apply Lemma [3~T1 Lemma l3~2l take the limit e \ and get 

m oo I , \ 

1=0 fc=o i=0 ^ ' 

f i+1 i oo ^ 

( T 2*c)"- fc - 2 ^ C log-''[T2*^] - -^-j- ^ d r ( T 2^r . 

2=0 ' r=0,r^n-fc-2 ' 
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We reorder the two infinite sums to one infinite sum, which can be done because 
of the structure of the coefficients dk , d r of the exponential integral that lets us 
compare the new coefficients to the coefficients of the exponential series, and 
get the expression (|4.14|) . □ 

Next, we need a similar expansion for the derivative (j)' u (u). 

Lemma 4.4. Fori = 0, lo E M\{0} ; n > 3 and fixed u > 0, the first u- derivative 
of 4>ui {u) satisfies the expansion 

n 

^(u) = -iue~ iuJU + h (oj, u)+ ^2 ( 2iuj T W(2i^) log k {-2iu)h ljk {u, u) , 

i>j+k=l 

(4.16) 

where the functions ho,hijk E C n (R) with respect to lo. 

Proof. In order to prove this, we use the fact that </>^(it) satisfies for u > an 
integral equation analog to (|2.10[) 

/>oo 

0l(u) = -iu)e~ iujiu - / cos(w(u - v))Vq(v)^ u (v) dv . 

J u 

We estimate cos(w(m — v)) and its lo— derivatives for real lo and v > u > by 

|^Jcob(w(u-i;))| < (v-u) n < (2v) n , n E N . (4.17) 

Thus, using this estimate for n = together with the iteration scheme (|3.14|) 
for ^) u] (u) 1 we obtain a well defined iteration scheme for the u-derivative: 

oo 

</>L( u ) = where 

fc=0 

V&°>(u) = -iue-^ u = UW)'(u), (4.18) 



^ +1 \u) = - cos{lo(u - v))V {v)^{v) dv = (^ +1 >) («) , 



with fc > 0. Due to this iteration scheme together with the estimates (|4.17p 
that replace the bounds (|3~T0| and the identity cos(w(u - v)) = 1/2 (e"^"^ + 
e iuj(«— u)^ we can use ^.^g decompositions of the , which we have made in 
the proof of Lemma I4TT1 In particular, we apply the procedure of this proof, in 
order to show the claim. □ 

Now, we use the expansions (|4.2l) . (|4.16p . in order to analyze the w-depen- 
dence of the essential part of the integral kernel 



Im 



4>u{u)(l>u{v) 



At this stage, it is enough to set n = 4 in (|4.2j) . (|4.16[) for our purposes. Look- 
ing at the integral representation (|2.5[) of the solution, we see that u E M is 
fixed while v E K varies in a compact set, the support of our initial data ^Sq. 
Due to the Picard-Lindclof theorem and the analytical dependence in lo of the 
Schrodinger equation from the coefficients, the expansions (|4.2|) . (|4.16j) extend to 
any u, and v, respectively, on compact sets. Moreover, the following properties 
follow directly by the construction of the expansions. 
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Corollary 4.5. For 4>i=j+k>l the function gijk, hijk can be constructed 
such that they obey the equalities 



g i]k (uj,u) + o(lo k ) = djk (e %uu + g (u,u)) and 

fnr i pvpti 

(4.19) 



hijk(u, u) + o(uo K ) = Cijk h (u;, u) , for i even 



g i3 k(u,u) + o(lo k ) = Cijk {e lu}u + g (u>,u)) and 
hijk(u>,u) + o(u> K ) = c ijk h (uj,u) , for i odd. 

where k is an arbitrary integer and the c^k are real constants, in particular not 
depending on u. 

Proof. We show this exemplarily for the first terms gno, ^no- hi this situation, 
(|4.19p holds because the first term, where (2iui) log(2iu;) appears, appears with 
C30/2 e luJU and there are no other terms with this w-dependence except the terms 
that are generated by this [cf. the calculations (I4.4|) . (I4.10|) ]. Thus, guo(uJ, u) 
is generated by e' luJU , which is just the complex conjugate of e -4 "", and this 
behavior is kept by the iteration scheme. So any C 4 -term that is generated 
is the complex conjugate of a corresponding term of go- This is valid, until 
one finishes the iteration scheme with the arguments at the end of the proof of 
Lemma 14. 1\ by what the o(w K )-term arises. Since one can do arbitrary many 
calculations and in each iteration at least a ±2iwlog(±2ia->) is generated, the 
K can be chosen arbitrary. Moreover, looking at the iteration scheme (|4.18p . 
the equalities for hno(uj,u) are a consequence of the arguments for gno(w, u), 
because of the fact that by the calculations concerning this scheme no additional 
highest order log-terms, i.e. i — j + k, are generated. □ 

In the following assume that k = 5. We expand the functions gijk(w, u) and 
hijk(&, u) in their Taylor polynom with respect to to at oj = up to the fourth 
order: 

4 i 

g ljk (u,u) = 2J — ] d ^9ijk\ {0 ^ m +r ijk (uJ,u) , 



m 

^ 1 



h ljk {u,u) = 2^ — d™h ijk \ (0u) LJ m + q ijk (uj,u) 



m 

m=0 



where the remaining terms rijk(w,u),qijk(u>,u) € C (R) behave for small u> as 
o(\lj\ ). Note that, due to this fact, any logarithmic irregularity multiplied with 
fijk,Qijk yields a C 4 -term with respect to u>. Moreover, we expand for fixed u 
the fundamental solution 4>u,(u) and its u— derivative 4>' u {u) 



k=0 k=0 

which exist, because these are analytic in u> for fixed u. Since the fundamental 
solutions q>,<j> are real for u> = 0, the coefficients <7o(0, u), ho(0, u), Cq(u) and 
do(u) are real for all u£l. Using all these properties, we expand 

w(<j>,<f>) 
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with the ansatz of a geometrical series with respect to ui. Note that, according 
to a result in [13 Section 6] the Wronskian does not vanish for uj — 0. By a 
straightforward calculation it is shown that, essentially using (|4.19p . the terms 
with the highest logarithmic order, i.e. (2iui) 1 log(2iujy \og(~2iuj) k ,i = j + k, 
vanish. Thus, we have to pick out the terms (2iu>) 2 log(2zw) J 1 log(— 2iuj) k with 
j + k = 1, in order to get the lowest regularity. Looking at the calculations (|4.4|) 
and (|4.10p [Note that these are the only possible terms, where a term with this 
irregularity appears the first time, according to our construction. The others 
are just a consequence out of these and hence a contribution to functions fltojjfe], 
the desired terms appear in the first time as 

e~ iuJU ((2iuj) 2 log(2iLu)(c + clogu) - c(2iuj) 2 \og(-2iuj)) , (4.21) 

where a (2ioj) 2 \og(2ioj)logu shows up in the first line of (|4.10|) . if one sep- 
arates log(— 2iu>u) = log(— 2iu>) + log it. All other such terms appearing in 
the second line of (|4.10p as well as in the second line of (14. 4|) vanish because 
of their coefficients. Applying the same arguments as before, it follows that 
#20i(w,w) + o(lu k ) = c(e~ tuu + g (uj,u)) and h 20 i(uj,u) + o(uj k ) = ch Q (uj,u). 
Hence, the terms with (2iui) 2 log(— 2iuS) cancel in the w-expansion of (|4.20p . 
Because of the additional log it-term, we get 

3210 (w, u) + o(uj k ) = ci{e-^ u + g (uJ, u)) + c 2 (e~ iuJU log u + g(u, «)) , 

and 

h 2W (uj, u) + o{u> K ) = cih (u, u) + c 2 h(uj, u) + -c 30 e lwu , 

with appropriate real constants c\,c 2 , where the last term appears by a di- 
rect calculation of ip^{u) with the part c^o/v 3 of the potential Vq(v). Fur- 
thermore, g(u,u),h(oJ,u) are C 4 -functions with respect to ui, where g(cu,u) is 
generated by the iteration of e _i ""logM and h(ui,u) the consequence out of 
this in (|4.18p . One directly verifies that g(0, u), h(0, u) are real, in general non- 
vanishing. Putting all these informations together, one sees that there appears a 
term with (2iui) 2 log(2iu;) in the w-expansion of (|4.20|) . which is generated on the 
one hand by the g(0, u), h(0,u), and on the other hand by the 2iu> log(2iw)-part 
multiplied with the w-contribution of first order of <j>,4>'. This represents the 
part with the highest irregularity with respect to uj. Moreover, the related co- 
efficients are purely real, depending on u, v and in general non-vanishing. Using 
the identity 

7T 

log(2iw) = i- sign(w) + log(2|w|) , 

and taking the imaginary part of (14.20p . which is just the essential part of our 
integral kernel, we obtain as the lowest regular w-term in the expansion of (|4.ip 
at uj — 

co{u)g 2 o{v) uj 2 sign(w) , (4.22) 

where the function g 2 o(v) arises out of the foregoing calculation. The symme- 
try of (|4.ip with respect to u,v yields immediately g 2 o(v) — kco(v) with an 
appropriate constant k ^ 0. 

In the next step we want to use (|4.22|) , in order to derive the decay of the 
solution *i>(t, u) given by (|2.5p . To this end, first we have to analyze the behavior 
of the w-derivatives of the integrand up to the fourth order for large \uj\. 
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Lemma 4.6. For u G M and compactly supported smooth initial data G 
Cq^R) 2 of the Cauchy problem, the u -derivatives of the integrand in the integral 
representation i2. 5}) 



d': 



Im 



,(u)<j» u (y) \ ( to 1 



'supp* \ w(9uiM j 

have arbitrary polynomial decay in lu for \u>\ — > oo. 



Vo(v)dvj , me {0,...,4} , 

(4.23) 



Proof. We proceed essentially as in the proof of [HI Theorem 6.5]. To this end, 
we have to investigate the behavior of <^ w (ii), 4>u>(v) in ui for u G R fixed and v 
in the compact set supp'I'o- We start with (j)^. We assume that \lo\ > 1 and 
Mo G R is arbitrary. Obviously we find for any v > u > u and m G {0, 4} a 
constant Ci(uo) such that 



3': 



1 



sin(w(w — w)) 



< J- Cl K)(i + H)' 



(4.24) 



Furthermore, splitting the potential as 

5 p-3 



\ — • \ — • log 9 ^) 



p=3 9=0 



and following an analog calculation as in (|4.6p . we obtain for the w-derivatives 
of the first iteration <j)u\u) for all u > 1 and m G {0, ...,4} the estimate 



< —c 2 u m - 2 



(4.25) 



with an appropriate constant C2. [Note that this is just an analogue to the 
estimate (|4.7p .] For all u < 1 and m £ {0, ...,4} we get 



< 



d" 



1 1 



+ 



00 y 

id 



sm(uj(u - v))V (v)e- luJV dv 

sin(oj(u- v))V {v)e- MV dv 
I 1 - 



where / is a continuous function with respect to u and the second term arises 
by the same method as we used for the estimate (|4.25p . Defining C4 by 



Ca 



max max j f/(m, u) + C 3 ±- jr \u\ k ) (l + M) 2 -" 1 ! 
me{o,...,4}«e[«o,i] [ \ M£^ J J 



and C5 := max(C2, C4), we obtain for all u > uq and m G {0, 4} the bound 
d^J\u) <-L C r (l + \u\r- 2 . (4-26) 
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In order to estimate the derivatives of the second iteration tj>}j (u) up to the 
fourth order, we subtract the first exact term out of the integration by parts in 

em 



C30 



Aioju 2 
m < 4 the bounds 



e luu , from the first iteration <pjj\u) and obtain for u > 1 and 



C 30 ^—iuu 



m— 3 



4tum 2 

Thus, in order to estimate the aj-derivatives of the second iteration: 

c 30 —iuv 



f 1 / 



+ 



d" 



AlLUU 2 



- wa.(u(u - v ))V {v)-^-^e- tuJV dv 



(4.27) 



dv 



Using the estimates (|4.27l) . (|4.24[) . and once again the method of splitting up the 
potential and integrating by parts for the second integral, we get for u > 1 and 
rn < 4 the bounds 

1 



\uj 



m — 4 



and thus, following the foregoing arguments for all u > uq (after possibly en- 
larging C5) the estimates 



d™€\u)) <-Lc 5 (i + \u\y 



(4.28) 



Using (j424]) and (|3~28"]) . we obtain for the ^-derivatives of the third iteration 
for all u > uq 

9^(u)\ < |E(T) jf d ™~ k (^M^-^)))v (v)d^Hv)dv 

< lediuo)^-- / (1 + \v\) m - 4 — V (v) dv . (4.29) 

Note that interchanging the integral and the w-derviatives is permitted, be- 
cause the w-derivatives of the integrand are integrable due to the estimates 
(|4.24|) , (|4.28p and the l/u 3 -decay of the potential V (v). We show by induction 
in n for all u > uq the inequality 

< 16CiS)C 5 pQ M (m, U ) ^ P u «- 3 , Vn>3, (4.30) 
where the functions Q UJ (m, u) and Pu,(u) are given by the integrals 



Q u (m,u) := 

P u (u) := 16d(u )C 6 



(l + \v\) m - 4 — V (v) dv 
\u>\ 



1 



V (v) dv 



where Cq is a constant chosen such that for all x > v > uq 

(1 + |a;|) fe - m < C 6 (l + \v\) k ~ m , < k < m < 4 
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The initial step is now given by (|4.29p . So assume that (|4.30p holds for n. Then, 
according to the iteration scheme, 



oL n+1) («) 



< 



OO -i 

Cl (u )(l + \v\r- k — Vo(v) 

\U!\ 



fc=0 

xlZC x {u )C 5 ±-Q u {k,v) - 1 

M (n-3)! 



Using the inequality 

Qu(k,v)<C 6 (l + \v\) k - m Q UJ (m,v) 

and the monotonicity of Q^, we obtain 

1 f°° 1 
d^ +1) (u) < 16Ci(«o)C 51 -r / WC^il + lvir^ — Voiv) 



xC 6 (l + \v\) k - m QUm,v) 



(n-3)! 

1 f°° dP 1 

< \QC 1 {u a )C 5 —M m ,u) / -t^(») 7 —PUv) n ~ 3 dv 

M y u aw (n-3)! 



16Cl(tto)C 5 ,— rQw (W, ti) ^TTT-Pw (w) 

|w| (n-2)! 



n-2 



and (|4.30jl follows. In particular, we get for all u > uq and m < 4 the estimate 



CiM-s;^ u " < c 5 ^(i + | u |) m - 2 + ^c 5 (i + | u |r- 4 



+16Ci(mo)C5t- -,Qui{m, u)e 



(4.31) 



and the right hand side obviously tends to zero as \u>\ — > oo. 

In an analog way using the iteration scheme (|2.12p for (j>, one shows for all 
u < uq and m € {0, ...,4} 

d™j>„{u) - 9™eH < Y -Mm, u) n = e M "(™>") - 1 , (4.32) 

I * — ' TV. 

n—1 

where M u (m, u) is given by 

M u (m,u) (1 + \v\) m V (v) dv , 

with a sufficiently large constant Cj. Note that this integral is well defined, and 
in particular the estimate is obtained easier, due to the fact that Vo(v) decays 
exponentially as v — > — oo. Moreover, the right hand side in (|4.32[) also goes 
to zero as \u>\ — ► oo. Thus, due to (|4.31j) and (|4.32p . the w-derivatives of the 
fundamental solutions up to the fourth order c*™<^(u),<9™0 w (?;) are controlled 
for large \lu\ by constants, which depend on u and the support of the initial data 
^o- One also shows with these results and applying the same arguments to 



that the Wronskian 



behaves as 0(M) and 



, m < 4 is 
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bounded by constants as \u>\ — > oo. Hence, interchanging in the representation 
(|4.23p the differentiation with respect to u> and the integral, which is no problem 
because of the compact support of ^q, making the substitutions 

dMv) = ^ + V (v)} u (v)) + ^ (-dj'^v) + V (v)dJ u (v)) 

as well as the analog substitutions for the second, third and fourth w-derivative 
[Note that in the region \ui\ > 1 4>uj(v) is C 4 with respect to u>, cf. Lemma T4.ll 
n = 4] and integrating by parts with respect to i>, one immediately has decay 
at least of 1/lo 2 . Thus iterating this procedure, which can be done because Vq 
and are smooth, yields arbitrary decay in lo and the lemma is proven. □ 

Remark 4.7. Since the method of the proof does not depend on the highest 
order w-derivative, the statement of Lemma 14.61 can be extended to arbitrary 
to. The only point where one has to be careful is the derivation of (|4.28p . since 
for ^-derivatives of higher order one has to calculate and subtract more exact 
terms than in (|4.27[) , due to convergence problems. If (14.28)) is not sufficient, 
in order to start the induction, one has to iterate this procedure appropriately 
many times. 

We are now ready to state and prove our main theorem: 

Theorem 4.8. Consider the Cauchy problem of the scalar wave equation in the 
Schwarzschild geometry 

□0 = 0, ((t>o,id t 4>o)(0,r,x) = $ (r,x) 

for smooth spherical symmetric initial data $o S Cq°((2M, oo) x S 2 ) 2 which 
is compactly supported outside the event horizon. Let = ((j)(t) , idt<j>{t)) € 
C°°(M. x (2 A/, oo) x S 2 ) 2 be the unique global solution which is compactly sup- 
ported for all times t. Then for fixed r there is a constant c = c(r, $o) such that 
for large t 

W)\ < I ■ (4-33) 

Moreover, if we have initially momentarily static initial data, i.e. dt4>a = 0, the 
solution <f)(t) satisfies 

W)\ < £ ■ (4-34) 

Proof. First, we decompose our initial data $o into spherical harmonics. Due to 
the spherical symmetry we obtain $o(?", ^, = $o( r )^oo(^j ( / 3 ) J where $o( r ) € 
Cq° ((2Af, oo)) 2 . Introducing the Regge- Wheeler coordinate u{r) and making 
the substitution ty(t,u) — r(u)<&(t,r(u)), our solution has the representation 

r 

where ^>(t,u) satisfies 
*(t,«) = 

- 1 ^ { I 4 ^#¥t) f £ I ) d -< (4 - 35) 

I" JR \7supp* \ W{(j)u,<t>u>) J V U U J J 



2G 



with initial data 4 , o(m) := r(u)$o(u) and the Jost solutions </>,</> in the case 
I = 0. According to the detailed analysis of (|4.20[) with respect to ui, the term 

Im( ^4r^T^T I ~ c o(«)52o(w)w 2 sign(w) - c 32 {u)g 32 (v)uj 3 log 2 |w| 



-c 3 i{u)g 3 i(v)uj 3 log \lu\ ~ c 30 (u)g 30 (v)u) 3 sign(w) 



is C 3 (K) with respect to w for fixed u G K, v G supp^o, where the Cij{u), 
gij(v) denote the appropriate coefficient functions. [Note that these are linearly 
dependent due to the symmetry of (|4.ip with respect to u,v.] Thus, defining 

/(«,«) := ( / **( f W 2 M * (^ , 

where the subscript denotes the first vector component, the term 

f(u,u) :=f(u),u)- ^c o (w)c?2o('0 2 )w 2 sign(w) +032(^)^32(^0) ^ 3 l°g 2 M 
+ c 3 1 (u ) d 3 1 (4>l ) u 3 log I lo I + c 30 (u ) rf 30 (-00 ) sign ( u> ) J 77 ( w ) 

is also C 3 (R) with respect to w. Here, ip^ denotes the second component of the 
initial data 'J'qj 



dij(A) : = / 9i](v)ip {v) dv , 

J supp ^/o 

and rj(u>) G Cq°(R) is a smooth cutoff-function which is identically to 1 on a 
neighborhood of ui = and outside a compact set. Moreover, because of 
Lemma |4~B1 the d™f(uj, u), m G {0, 1, 2, 3} have rapid decay for large \u>\ and are 
in particular L 1 (R) with respect to u. Thus, due to (|4.35|) . the first component 
of \& satisfies 

i> l (t,u) = — [ e- wt f(u,u)duj- - [ e- tut r{oj,u)du 
7T Jr n Jr 

I 



{it) 3 ir 



f(uj,u)d 3 e- luJt dio+ / r(u>,u)d 3 e- lult dw 



We write the second integral as /_ + Jq°°, integrate every integral three times 
by parts and obtain 

i>\t,u) = — ^- (ic (u)d 20 (^)+ [ e- iut d 3 J(tJ,u)cLo 
W d 7r V Jr 

pQ poo 

+ / e- wt d 3 r(uj,u)dtu+ / e' wt d 3 r(uj, u) duj 



Note that the other boundary terms vanish, because the d™ f(uj),m < 3 have 
rapid decay and t](ui) = outside of a compact set. Obviously, all integrals are 
well defined, and the Riemann-Lebesgue lemma shows the claim in the first case. 
If the initial data is initially momentarily static, all the di 3 (0 2 ) vanish and the 
entries in the matrix in (|4.35|) yield an additional uj. Hence, the highest irregular 
term is co(u)d2o(V'o) u;3 S1 g n ( a -')i an d the same arguments as before conclude the 
proof. □ 
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Remark 4.9. The decay rates 1/i 3 , and 1/i 4 , respectively, are optimal in the 
sense that there exists initial data such that these cannot be improved. This is 
obvious due to the fact that cq(u) > 0. 



5 Discussion on the case / ^ 

According to Price's Law [IS], the Zm-component & m (t, u) = ^W m (t,u) of a 
solution for the Cauchy problem for the scalar wave equation in Schwarzschild 
spacetime with compactly supported smooth initital data generally falls off at 
late times t as t~ 2l ~~ z and t~ 2l ~ 4 for initially momentarily static initial data, 
respectively. This has been confirmed in the previous section for spherical sym- 
metric initial data, i.e. in the case I = [cf. Theorem 14. 8j . Moreover, there is 
numerical evidence which lets us conjecture this to be correct [H]. We briefly 
discuss whether the methods of the preceding section still apply to the case 
when the angular mode I is non-zero. 

To this end, let us reconsider the construction of the fundamental solutions 
<j)ui °f the Schrodinger equation (|2.6p . First, we make some remarks about 
the fundamental solutions u> <j) u i{u) (see also [TSJ Section 5]). The fundamental 
solutions were constructed as the series 



lj%(u) = 5>£ m) N, (5.i) 

m— 

where the (j)^ are given by the iteration scheme 

cf>i m+1 \u) = - / S u (u,v)Wi(v)^(v) dv , (5.2) 



with potential, cf. also Lemma 



Wi(u) = Vi(u) — = c 3 i — - — I T + h{u), (5.3) 



where h{u) = O 



log 2 u 



for large u, and Green's function 

Su(u,v) = (h 1 (l,Ldv)fi2(l 7 uJu) — hi(l,Ldu)h 2 (l,uJv)) , (5.4) 



where 



hi{l,uu) = ^^J 1+1/2 (uju) , h 2 (l,uu) = y^pJ_j_i/ 2 (wu) 



(5.5) 



and J v denotes the Bessel function of the first kind. As initial function ififfl (u) 
we have chosen 

^) { u)=Je-^m^ H % 1/2 {uu), 

(2) 

where H„ ' denotes the second Hankel function. Since I is an integer, these 
functions are directly connected to the spherical Bessel functions and simplify 
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significantly. Namely, hi , hi have the following representations [cf. [TJ Chapter 
10]] 

hi (I, uju) = P(l + — , uu) so\(uju — + QQ + 2 ' UJU ^ cos ( UJU ~ ^' 7r ) (^-^) 
wm) = P(/ + — , uu) cob(uju + 2^ 7r ) — 2' sin(o;u + iz7r) (5.7) 

where P, Q are finite polynomials given by 

m 



fc=0 



with 



2' y fc!r(;-fc + i) 

And the initial function can be expressed by 

Due to the recurrence formulas for the derivatives of the Bessel functions, we 
have the identities 

d u h\{l,uju) = uh\(l — l,um) hi(l,u>u) , 

UJ 

duhztiiWu) = —uhzil — l,ijju) h<2(l,u)u) . 

UJ 

As a consequence, 

d u S u (u,v) = - ^ + 1 S u 

UJ 

+ v (hi (I — 1, Luv)h2(l, uju) + /ii(/, u)U)h2(l — 1, 

(-1)' 

+ u (hi(l,u>v)h2\l — I, urn) + hi(l — l,oju)h2(l,uiv)) . 

UJ 

This allows us to derive the necessary estimates for the Green's function S w (u, v). 
Exploiting the asymptotics we have already seen in Section 5] 



\suu,v)\<ci(—^--) (j-^rj-) 

\l + \uj\uj \l + \uj\vj 



l+l 

v\ Im aj|+nlm uj 



for v > u > and an appropriate constant C\. In order to derive an estimate 
for duSu and small \uj\, we make use of 

hi(l,uju) ~ ki(uju) l+1 + k 2 (uju) l+3 
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h 2 {l,Luu) ~ k 3 (ivu) 1 + fc 4 (um) 1+2 , if|w|u<Cl, 

and certain constants hi, ki [refer to the series expansion of the Bessel func- 
tions [TJ 9.1.10]] to obtain (note that v > u > 0), 

For large arguments \lu\u 3> 1 we use (|5.6j) . (|5.7|) and get by a straightforward 
calculation 



—21 

duiS u (u,v) ~ — 5- sin(w(u - v)) + d u 



Together with (|3.10|) . we obtain 



— sin(ui(u — vj) 

LU 



if \oj\u > 1. 



\S u (u,v)\ <C 3 - V ] , e "' Im "l +uIm " , if HOI. 
Combining these estimates, we find a constant C such that 

IWu,«)|<c(— ^ (V- ^j-pV »e"l to " , (5.9) 
\1 + / \1 + y 



for d > u > 0. Moreover, looking at (15. 8|) we get the following bounds for the 
initial function, 



|0i o) («)| < C 4 f— H-p) e" Im ", (5.10) 

|^2V)I < C 5 (—^—) ue^. (5.11) 

These estimates allow us to proceed in exactly the same way as in the proof of 
Lemma 13.41 As analogon to (u) we obtain the term 



which we calculate using (|5.6|) , (|5.7[) and (|5.8p . Essentially, we get integrals of 
the shape 

(„ ^ r ^ ^ v ^ dv + ^ jogv dv 



(wm)"!!)" 1 ^^ 1 V J v 3+k+m J v 3+k+m 

where q £ {0, 1}, < n, m, k < I. Note that the terms involving lu singularities 
resolve, due to the fact that u q> u is continuous with respect to u>. Computing 
these integrals via Lemma [3~T1 (in the limit e — > 0), we see (as before) that the 
only terms not being C 1 with respect to u> are of the form 

e lulu ujm l +k+1 {2iu J ) k+m+2 (log 2 (2zc^) + logulogpiwu) + log(2»w«)) , (5.12) 

modulo coefficients. Now, we apply the same iteration with analog estimates 
and all in all we have shown: 
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Lemma 5.1. For I > 1, u> 6 K\ {0} and fixed u > i/ie fundamental solutions 



where the functions 53,94,55 and 96 « r e C 1 (R) with respect to uj. 

Hence, we still have finite expressions for the Green's function S u (u,v) as 
well as for the initial function 4>^(v), which involve essentially the plane waves 
e w ° , e IWV . Expanding all these expressions and deriving estimates analog to 
(|5.9p and (|5.11jl for higher order w-derivatives, we can improve Lemma 15.11 in 
the same way as Lemma 13.41 following the arguments of the proof of Lemma 
14.11 Also, a similar result to Corollary 14. 51 seems straightforward. The problem 
now arises, when we have to derive an uj— expansion of the essential part of the 
integral kernel 



The main difficulty can be seen as follows. If we proceeded in the same way as 
in the case 1 = 0, the lowest regular term with respect to uj should appear with 
the power uj 21+2 [cf. proof of Theorem 14.8] in order to satisfy Price's law. But 
due to the fact that the first irregularity in uj looks as follows, 

e luJu u~ l 2iuj(c\og 2 (2iuju) + c\ogu\og{2iuju) + c\og{2iujuj) , 

[cf. equation (|5.12| ], we would have to find a systematic way in order to check 
that the coefficients in front of terms with lower regularity vanish. Because of the 
complexity of the calculations we did not succeed in this point. Thus, following 
the same arguments as for I = together with the analog result to Corollary 
I4.5[ which would involve 2iuj log 2 (2iuj) as highest irregularity, we would have to 
assume uj log |w| as the lowest regular term in the expansion of ()5.14j) . Except for 
this problem, we do not expect any further difficulties in extending Lemma 14.61 
to I 7^ 0, apart from the complexity of the calculations and the estimates. Thus, 
for arbitrary I it follows a similar statement to Theorem l4.81 but with the decay 
< c/t 2 , and in the case of momentarily static initial data \4>(t)\ < c/t 3 , 
respectively. The proof uses essentially the arguments of the proof of Theorem 
14.81 with the difference that one basically has to check the inequality 



To this end, one makes the substitution z — ujt and splits up the integrals to 
obtain 




(5.13) 




(5.14) 
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Computing the second integral and integrating the last two integrals by parts 
yields 

= \(J \og\z\e~ tz dz + ^ J - z e~ lz dz + 1 jf ^-e~ lz dz^j , 

and the inequality follows, after having integrated the last two integrals once 
again by parts followed by standard integral estimates. However, in view of 
Price's law, this result is not satisfying. 

Acknowledgments: The author would like to thank Felix Finster, University 
of Regensburg, for introducing to this problem and also for his helpful sugges- 
tions and remarks. 
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